A Successive Linearization Method Approach to Solve Lane-Emden Type of Equations

We propose a new application of the successive linearization method for solving singular initial and boundary value problems of Lane-Emden type. To demonstrate the reliability of the proposed method, a comparison is made with results from existing methods in the literature and with exact analytical solutions. It was found that the method is easy to implement, yields accurate results, and performs better than some numerical methods

Computing Simple Roots by an Optimal Sixteenth-Order Class

The problem considered in this paper is to approximate the simple zeros of the function by iterative processes. An optimal 16th order class is constructed. The class is built by considering any of the optimal three-step derivative-involved methods in the first three steps of a four-step cycle in which the first derivative of the function at the fourth step is estimated by a combination of already known values. Per iteration, each method of the class reaches the efficiency index , by carrying out four evaluations of the function and one evaluation of the first derivative. The error equation for one technique of the class is furnished analytically. Some methods of the class are tested by challenging the existing high-order methods. The interval Newton's method is given as a tool for extracting enough accurate initial approximations to start such high-order methods. The obtained numerical results show that the derived methods are accurate and efficient

The Mixed Finite Element Multigrid Method for Stokes Equations

The stable finite element discretization of the Stokes problem produces a symmetric indefinite system of linear algebraic equations. A variety of iterative solvers have been proposed for such systems in an attempt to construct efficient, fast, and robust solution techniques. This paper investigates one of such iterative solvers, the geometric multigrid solver, to find the approximate solution of the indefinite systems. The main ingredient of the multigrid method is the choice of an appropriate smoothing strategy. This study considers the application of different smoothers and compares their effects in the overall performance of the multigrid solver. We study the multigrid method with the following smoothers: distributed Gauss Seidel, inexact Uzawa, preconditioned MINRES, and Braess-Sarazin type smoothers. A comparative study of the smoothers shows that the Braess-Sarazin smoothers enhance good performance of the multigrid method. We study the problem in a two-dimensional domain using stable Hood-Taylor Q2-Q1 pair of finite rectangular elements. We also give the main theoretical convergence results. We present the numerical results to demonstrate the efficiency and robustness of the multigrid method and confirm the theoretical results

Thermal Radiation and Buoyancy Effects on Heat and Mass Transfer over a Semi-Infinite Stretching Surface with Suction and Blowing

This study sought to investigate thermal radiation and buoyancy effects on heat and mass transfer over a semi-infinite stretching surface with suction and blowing. Appropriate transformations were employed to transform the governing differential equations to nonsimilar form. The transformed equations were solved numerically by an efficient implicit, iterative finite-difference scheme. A parametric study illustrating the influence of wall suction or injection, radiation, Schmidt number and Grashof number on the fluid velocity, temperature and concentration is conducted. We conclude from the study that the flow is appreciably influenced by thermal radiation, Schmidt number, as well as fluid injection or suction

Numerical Solution of Nonlinear Fractional Volterra Integro-Differential Equations via Bernoulli Polynomials

This paper presents a computational approach for solving a class of nonlinear Volterra integro-differential equations of fractional order which is based on the Bernoulli polynomials approximation. Our method consists of reducing the main problems to the solution of algebraic equations systems by expanding the required approximate solutions as the linear combination of the Bernoulli polynomials. Several examples are given and the numerical results are shown to demonstrate the efficiency of the proposed method

Convergence Analysis of Legendre Pseudospectral Scheme for Solving Nonlinear Systems of Volterra Integral Equations

We are concerned with the extension of a Legendre spectral method to the numerical solution of nonlinear systems of Volterra integral equations of the second kind. It is proved theoretically that the proposed method converges exponentially provided that the solution is sufficiently smooth. Also, three biological systems which are known as the systems of Lotka-Volterra equations are approximately solved by the presented method. Numerical results confirm the theoretical prediction of the exponential rate of convergence

Two Optimal Eighth-Order Derivative-Free Classes of Iterative Methods

Optimization problems defined by (objective) functions for which derivatives are unavailable or available at an expensive cost are emerging in computational science. Due to this, the main aim of this paper is to attain as high as possible of local convergence order by using fixed number of (functional) evaluations to find efficient solvers for one-variable nonlinear equations, while the procedure to achieve this goal is totally free from derivative. To this end, we consider the fourth-order uniparametric family of Kung and Traub to suggest and demonstrate two classes of three-step derivative-free methods using only four pieces of information per full iteration to reach the optimal order eight and the optimal efficiency index 1.682. Moreover, a large number of numerical tests are considered to confirm the applicability and efficiency of the produced methods from the new classes